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Abstract. Recently obtained recurrence formulae for relativistic 
hydrogenic radial matrix elements are cast in a simpler and perhaps 
more useful form. This is achieved with the help of a new relation 
between the r a and the (3r b terms (j3 is a 4 x 4 Dirac matrix and a, b are 
constants) in the atomic matrix elements. 
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The relation of quantum calculations with experimental results is usually made at 
the level of expectation values or matrix elements of appropriate operators (Moss 1972, 
Wong and Yeh 1983a,b, De Lange and Raab 1991, Quiney et al 1997). In atomic or 
molecular physics this means matrix elements of powers of a radial coordinate between 
hydrogenic states. This happens both because the hydrogen atom is exactly soluble 
— hence very useful in approximate calculation schemes — and because a power of r 
can be either a well defined function of interest, like in the London or the Lennard- 
Jones potentials, or as a term in the multipolar expansion of the electromagnetic 
field. Given such importance, we have been working in obtaining recurrence relations 
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between matrix elements of radial functions between hydrogenic eigenstates but in 
the relativistic realm (Martinez-y- Romero et al 2000, 20001). Our purpose here is 
to report on a new relationship we have recently found between hydrogenic matrix 
elements of powers of r and of powers of j3r that can simplify the task of using the 
recurrence relation previously reported. 

By extending to Dirac's relativistic quantum mechanics a hypervirial method 
(Nunez- Yepez et al 1995) that is known to led to the nonrelativistic Blanchard 
recurrence relation (Blanchard 1974), it has been possible to obtain the relativistic 
recurrence formulae (Martinez-y-Romero et al 2000, 20001) 
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b =4A [(£-) 2 "4c 4 ], 
b 2 =c 2 (l - A) [(A 2i r - 4A 2 ] , 
e =2E + [E-A 21 -4c 2 A], 
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(5) 



e 2 =c 2 ^[(A 21 ) 2 -4A 2 ]. 



Here we defined E ± = E 2 ± E±, f3 = diag(l, —1) is a 2 x 2 Dirac matrix, the kets 
|1) and 1 2) stand for radial eigenstates of the relativistic hydrogen atom, solutions 
of the Dirac equation H k ip k (r) = E k ipk(r) 7 k = 1,2, where k is a label that will 
serve to distinguish the states appearing in the recurrence relations, i/jk(r) is a radial 
bispinor, and A^ = e 2 (2j 2 + 1) ± ei(2j'i + 1). The quantum number e equals +1 
when / = j + 1/2 and equals —1 when I = j — 1/2 where j = 1/2, 3/2, 5/2, ... is the 
total angular momentum quantum number, and / is the orbital angular momentum 
quantum number. It is to be noted that the relativistic states are often labeled by the 
quantum number k — as done in (Drake 1996) — instead of by e; the relation between 
these two numbers is k = — e(j + 1/2). The possibly complex numbers A are the 
constant exponents of the radial coordinate r. In the previous equations as in all the 
paper we use atomic units (h = m e = e = 1). 

Written in full, the radial Dirac Hamiltonian Hk and the energy eigenvalues can 
be written as (Drake 1996, Martmez-y-Romero et al 2000, 20001) 



where p r = — (i/r)(l + rj^), a r is a 2 x 2 antidiagonal matrix having — l's as the only 
nonzero elements, (r) is a radial potential behaving as the fourth component of a 
4- vector. As it must be clear from (6), the function ipk(r) is 



where n k = 0, 1, 2, . . . is the principal quantum number, and = j\ ± 1/2, according 
to whether refers to the big or to the small component of the hydrogenic spinor, so 
Fk and Gk are, respectively, the big and the small component of the radial bispinor 
with quantum numbers e^, jk, and n^, corresponding to an energy (Greiner 1991) 
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where a F = 1/c ~ 1/137 is the fine structure constant. 

The recurrence relations (1) and (2) could be more useful if the matrix elements 
of r A were uncoupled from the matrix elements of f3r x . It is the purpose of this Letter 
to show how to achieve such uncoupling [see also the rather cursory discussion in 
(Ilarraza-Lomeli et al 2001)]. 

To evaluate the relation that is needed to uncouple the recurrence relations, let us 
first evaluate H 2 £ + £Hi, that is 
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where £(r) = H 2 f(r) — f(r)Hi, f(r) a radial function, and we indicate with primes 
r-derivatives. Next evaluate 
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taking (9) and (10) together, we get 
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we may get 
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We are almost done with the equation, but we need to eliminate the term in the left 
hand side of (13), to this end we have to evaluate the following 



H 2 f + fH l = -ica r 



2f^ + f' + -f + P^f 
ar r 2r 



+ 2f((3c 2 + V). (14) 



Letter to the Editor 



5 



The next steps are analogous to the previous ones, but now using equations (14) and 
(10); by further juggling with them, we get 
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Combining (13) with (15) and taking matrix elements between the states (s 2 \ and \si), 
we finally obtain the fundamental equation for the disentanglement: 
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where we have explicitly used both that f(r) = r x and the potential as Coulomb's 
V(r) = —Z/r. The coefficients are given explicitly as 



9i 



Ji 



,'A" A + 

2 / ^21^21 



k 21 



AJi + A+) 



E + + c 2 (X-l) 



(17) 



j 2 =-^Z(A 2 - 1 + A+). 

Equation (16) is a recurrence relation that can be found useful for various tasks 
but for now we just want to uncouple relations (1) and (2). 

We need first to equate r A taken from equation (2), to r A taken from equation (16); 
in the resulting equation there is yet a term involving r x ~ 2 which we eliminate using 
again (16). This gives us the following relation between matrix elements involving f3r x 
terms only, 



Vo(s2\(3r x \ Sl ) = Vi&lPr^ 1 ^) + 77 2 (s 2 |/3r A - 2 | Sl > + m (s 2 |/3r A - 3 | Sl >, (18) 

the explicit form of the coefficients is given after equation (19). 

To derive the recurrence relation involving matrix elements of r A terms only, we 
get the term [3r x ~ 2 that arises from equating the term /3r A_3 taken from (1) to the 
same term taken from (16). From here we also obtain also the terms f3r x and /3r A_1 
that are to be substituted into (2). Doing such substitutions we finally obtain the 
sought after recurrence relation 

^o(s2|r A |si) =z/i(s 2 |r A ' 1 |si) + z/ 2 (s 2 |r A ~ 2 |si> + ^(s2|r A - 3 |si> 

(19) 

+ ^4(s 2 |^ 4 |si) + v 5 (s 2 \r x - 5 \ Sl ). 
The coefficients in equations (18) and (19) are given by 
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where the symbols defined for writing the recurrence relations are: M = A(A — 1)E + , 
and 



L> =A-£;- - 4c 2 A, 
J(2M - A+D) 



J 



£> + 4c 2 



G = 
H = 
K = 



(X-1)D 
D 



D + 4c 2 ' 
S 



F = (E-) 2 -4c\ 
S = A++A 2 -, 



(22) 



^21 



L = 4A 2 - (A 



21/ ) 



The recurrence relations (16), (18) and (19) hold true if oji + uj 2 + 1 + |A| > 
(Martinez-y-Romero et al 2000, 2001) where the u a = + v / (j a + 1/2) 2 - Z 2 a 2 F , a = 
1, 2 are real numbers. This conditions is basically that any integrand goes to zero faster 
than 1/r as r — > oo. As it is easy to realize, the new recursions give no new information 
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in the case of the relativistic Pasternack-Sternheimer (1962) rule — i.e. Aj x = A = 0— 
a relation which must then remain as given in our previous contributions (Martinez- 
y-Romero et al 2000, 2001). The new uncoupled relations should be anyhow useful in 
relativistic atomic or molecular physics calculations. 

This work has been partially supported by PAPIIT-UNAM. We acknowledge useful 
discussions with A Ilarraza-Lomeli and N Valdes-Martinez, and the friendly support 
of Zura, Mec, Shapka, Rabita, Orejitas, Weiss, Gorbe, Bieli and all the gang. 
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